Non-thermal fixed points: effective weak-coupling for strongly correlated systems 

far from equilibrium 
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Strongly correlated systems far from equilibrium can exhibit scaling solutions with a dynamically 
generated weak coupling. We show this by investigating isolated systems described by relativistic 
quantum field theories for initial conditions leading to nonequilibrium instabilities, such as paramet- 
ric resonance or spinodal decomposition. The non-thermal fixed points prevent fast thermalization 
if classical-statistical fluctuations dominate over quantum fluctuations. We comment on the possible 
significance of these results for the heating of the early universe after inflation and the question of 
fast thermalization in heavy-ion collision experiments. 
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INTRODUCTION 

Important phenomena in early universe cosmology 
("Big Bang") and collison experiments of heavy nuclei 
(" Little Bangs" ) involve quantum fields far from equilib- 
rium. A prominent topical example concerns the role of 
nonequilibrium instabilities for the process of thermal- 
ization. The heating of the early universe after inflation 
may proceed via an instability such as parametric reso- 
nance 0, Q ■ Similarly, plasma instabilities may play an 
important role in our understanding of observations at 
the Relativistic Heavy Ion Collider [3(. Instabilities also 
arise in many other areas, such as dynamics of ultra-cold 
quantum gases. 

Nonequilibrium instabilities lead to exponential 
growth of field fluctuations on time scales much shorter 
than the asymptotic thermal equilibration time. Though 
their origin can be very different, the subsequent evo- 
lution after an instability follows similar patterns: Af- 
ter a fast initial period of exponential growth the dy- 
namics slows down considerably. At this stage all pro- 
cesses become of order unity and one is dealing with a 
strongly correlated system that has to be treated non- 
perturbatively, even if the underlying microscopic theory 
is weakly coupled. The subsequent evolution is charac- 
terized by power-law behavior reminiscent of turbulence. 
It has been argued that this behavior does not occur in 
the non-perturbative regime and a perturbative analysis 
is employed with apparent success [J, [I| . 

In this Letter we show that far-from-equilibrium dy- 
namics in the non-perturbative regime can approach scal- 
ing solutions with a dynamically generated weak cou- 
pling. As an example we consider scalar iV-component 
quantum field theory with quartic self-interaction fol- 
lowing a parametric resonance instability. In the non- 
perturbative regime we find new scaling solutions with 
strongly enhanced low-momentum fluctuations ~ 



as compared to a high-temperature distribution ~ p~ 1 . 
They correspond to non-thermal fixed points of the 
time evolution equations for correlation functions once 
classical-statistical fluctuations dominate over quantum 
fluctuations. At sufficiently high momenta we recover the 
perturbative behavior ~ p~ 3 / 2 reported in the literature. 

We employ the two-particle irreducible (2PI) large- N 
expansion to next-to- leading order The 2PI approx- 
imation schemes have been applied to a variety of far 
from equilibrium phenomena, including parametric reso- 
nance Q and tachyonic preheating Q. They are known 
to describe the late-time approach to thermal equilibrium 
characterized by Bose-Einstein or Fermi-Dirac distribu- 
tions, respectively i|. The non-perturbative regime af- 
ter an instability is traditionally described using classical- 
statistical field theory simulations, and we present a com- 
parison of quantum and classical evolution. 



NONEQUILIBRIUM INSTABILITIES 

We consider a relativistic real scalar iV-component 
quantum field $ a (a = 1, . . . , N) with A/(4!iV) (^A) 2 
interaction, where summation over repeated indices is 
implied. The macroscopic field of the quantum theory 
is given by 4> a {x) = (& a (x)), where the brackets describe 
the trace for given initial density matrix. There are two 
independent two-point correlation functions, which can 
be associated to the anti-commutator and the commuta- 
tor of two fields: 

F ab (x,y) = ^({*o(aO, *&(»)}) -Mx)Mv), 

p ab (x,y) = <[* a (aO,* 6 (i,)]>. (1) 

Here p is the spectral function, which is related to 
the retarded propagator Gn{x,y) — p(x,y)Q(xo — yo)- 
The statistical function F is proportional to "occupa- 
tion number", which may be taken as n(t, p) + 1/2 = 



2 



1/2 

[F(t,t';p)dtdt'F(t,t';p)] t L_ t , for spatially homogeneous 
systems [6(. 

The nonequilibrium time evolution of F and p is de- 
scribed by coupled differential equations, 



[□ x <5 ac + Ml c (x)] p cb (x, y) = / d d zY> p ac (x,z)p cb (z 1 y), 

Jxq 

rVo 

p x 6 ac + M* c (x)]F cb (x,y) = I d d z^ c (x,z)p cb (z,y) 



d d zZP c (x,z)F cb (z,y), 
(2) 

and a similar equation for (f> a (x) [6]. These evolution 
equations would be exact for known effective mass term 
M 2 (cj),F), spectral (imaginary) part of the self-energy, 
S p (0, F, p), and statistical (real) part, S F (0, F, p). Here 
ti describes the initial time and we will consider a pure 
initial quantum state with spatially homogeneous fields 



<t> a (t) =(T(t)y/6N/\6 al 



(3) 



for A <C 1. All quantities will be given in units of the 
initial ao = a(t = 0). The initial two-point functions are 
taken to be diagonal with F ab = diag{f|| , F±, . . . , F±} ~ 
O(N X°) and the initial spectral function is fixed by the 
equal-time commutation relations. 

Dynamics of parametric resonance in quantum field 
theory has been studied in detail Q using the 2P1 l/N 
expansion to NLO [6j. Primary resonant amplification 
occurs for the dominant transverse modes in the mo- 
mentum range < p 2 < a 2 /2. This is followed by a 
non-linear regime, where F± ~ 0(N 1 ^X~ 2 ^). Here oc- 
cupied low-momentum modes act as sources for a sec- 
ondary stage of enhanced amplification in a broad higher 
momentum range. The exponential growth stops when 
F± ~ Fit ~ O(N X^ 1 ) and the dynamics slows down 
considerably. At this non-perturbative stage all processes 
are of order unity and the system is strongly correlated. 

Fig. [T] shows Fx (t, t, p) for times t = 1 0, t — 40 in the 
non-linear regime, and t — 90 in the non-perturbative 
regime. The solid line shows the result for the quan- 
tum evolution for A = 0.01 and N = 4. For compari- 
son, the dashed line gives the same quantity as obtained 
from simulations of the corresponding classical-statistical 
field theory on a lattice with same initial conditions fol- 
lowing Ref. 0]. The precision of agreement between the 
curves is remarkable. Quantum fluctuations are expected 
to be suppressed if the classicality condition F 2 3> p 2 is 
fulfilled [9j. However, with p 2 being of order unity the 
quantum corrections turn out to be extremely small even 
for F 2 -c f • We emphasize that this is not a generic prop- 
erty of the NLO approximation 0, Q , but a consequence 
of the instability dynamics and A < 1. In turn, contri- 
butions beyond NLO seem to play an inferior role even 
for the non-perturbative regime. 
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FIG. 1: Two-point function F±(t, t;p) as a function of mo- 
mentum |p| for three different times. Quantum (solid) and 
classical (dashed) evolution agree well even where fj_(p) <C 1. 



NON-THERMAL FIXED POINTS 

Fixed points are time and space translation invariant 
solutions of the evolution equations (|2|), which require 



KAp)F cb {p) 



,(p) Pcb (p) = (J$(p) + ■/#(?)) /A 

(4) 



= 0. 



Thermal equilibrium solves ^ using the fluctuation- 
dissipation relation F^l q \p) = [nrfpo) + 1/2] P^ b \p) 
and correspondingly for self-energies [6| . We show Q has 
non-thermal approximate solutions if classical-statistical 
fluctuations dominate over quantum fluctuations. 

Firstly, we analytically determine fixed point solu- 
tions with F ab (p) = f(p)/XS ab , Pab(p) = p(p)S a b and 
a = const 7^ 0. Secondly, these non-trivial solutions are 
shown to describe well the slow dynamics in the non- 
perturbative regime. At NLO of the 2PI l/N expansion 
in the classical-statistical field theory limit we have 

4 3 a W) = ~5^4 jd A k^ q 8\p- q -k) 

[Aeff(fc) + A off (<z) + A cff (p)] [p(k)f(q)f(p) 
+f(k)p(q)f( P ) - f(k)f(q)p(p)} , (5) 
where we summed over components a. Here 

f 



XMp> n+iuwr 

with the resummed "one- loop" retarded self-energy 
I 



3(2tt) 4 



d 4 qf(q)G R (p-q). 



(6) 



(7) 



Diagrammatically, contains contributions described 
by the "one-loop" graph in Fig. [2J In contrast to per- 
turbative Feynman diagrams, here lines correspond to 
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FIG. 2: NLO contributions to the fixed point equation @. 

self-consistently dressed propagators F or p Q. The full 
dot denotes the effective vertex ((fjj), which is graphically 
displayed in the lower graph of Fig. [2] It iteratively gen- 
erates the infinite series of resummed "chain" diagrams 
contributing at NLO in the 1/N expansion. Similarly, 

J <£W) = Ys^r / d4fcd4 « d4 ^ 4 (f+ fc -9- r ) 

X cS (p + k){[f(p)p(k)+p(p)f(k)}f(q)f(r) 

-f(p)f(k)[f(q)p(r)+p(q)f(r)}} 
(l-A(p + k;a)) (8) 

corresponds to the "two-loop" graph in Fig. 2] with 
A(p; a) = 4a 2 Re[G fl (p)/(l + IL R (p))}. 

In the non-perturbative regime, where the dynamics 
is slow, one can extract quantities such as X c s(p) also 
directly from the nonequilibrium quantum evolution by 
Fourier transformation with respect to relative coordi- 
nates. The result is shown in Fig. [3] for t = 240 as a 
function of three-momentum p for different values of the 
frequency u) = po, with parameters as in Fig. [TJ For 
momenta larger than order unity X e g(p) tends to one. 
One observes that this holds if either p or w is not small. 
Therefore, in this range the dynamics is well approxi- 
mated by the one- and two-loop expressions as displayed 
in Fig. [2] without the effective vertex. Most strikingly, for 
small four- momentum the infinite series of 0(1) contribu- 
tions add up to a substantially reduced effective coupling. 
The latter is still slowly changing in time and is shown 
to approach a power-law behavior in the following. 

For scaling solutions the correlators behave under si- 
multaneous scaling p — > sp and <f> — > sip as [3] : 



F(p) 
P(P) 



„2+a 



F(sp) 
= s 2 p(sp) , 
= s 7 X c e(sp) . 



(9) 



Similarly, we have Gr(p) — s 2 Gn(sp) and it follows that 
= s a Tlji(sp). Taking sp ~ 1 one concludes for 
a > that U R (p) > 1 for p -> 0. According to © the 
infrared scaling behavior of the effective vertex in this 
case is described by © with the exponent 



7 = —2a . 



(10) 
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FIG. 3: The effective vertex at t = 240 as a function of spatial 
momentum p for different frequencies ui = pa- 



Using this scaling behavior we find 



s°jiV(sp;sa) and J a V(p;0) = s a JaV (sp; 0) such that 



Ja 3 J (p; o) = 



Jaa (p; 0) dominates the infrared dynamics, 
dependent part JaV (p; a) 
since A(p; a) = s~ a A(sp; sa) in this regime. 

Accordingly, for the infrared we consider solutions of 
/ d 3 pJaa (Po,P', 0) = 0. The integration over spatial mo- 
mentum allows us to solve this equation using generalized 
Zakharov transformations flfjj . In this way the problem 
is reduced to simple algebraic conditions for the expo- 
nents. E.g., to map the second term in the integrand 
sum of ((8]) onto the first we employ the transformation 



The field- 



JaV (p; 0) scales as J$ (p; a) 



P 



Po/ko, qo 



qoPo/ko, r — ► r a p /k as well as 
kpo/fco, k —> ppo/fco, q ->• qpo/fco, r — * rp /k , and 
similar for the third and fourth term. Using the scaling 
properties ([9]) the integrand vanishes if a = (4 — 7) /3 or 
a = (5 — 7)/3. With (TIT)]) this yields two non-thermal 
fixed point solutions for the infrared scaling behavior: 



a 



5. 



(11) 



We do not list classical thermal solutions (a = 1, 
0), which appear as a consequence of the fluctuation- 
dissipation relation discussed above. We emphasize that 
(jTTJ) cannot be obtained from a perturbative 2 <-> 2 scat- 
tering analysis 0, E3] , which we recover for vanishing 7. 

The size of the one-loop retarded self-energy ([7]) dis- 
tinguishes the low-momentum from a high-momentum 
regime. The latter is defined by Hr(p) <C 1 and 
in this range the effective vertex scales trivially (7 = 
0) according to ((6]) with \ e ff(p) — 1. Therefore, 

= s 2a JaV (sp; sa) and 



at high momenta JaV (p; a 



J i a i a > (p;0) 



a 3a t(4) 



Jaa(sp;0) such that JaV(p;a) domi- 



nates over JaV(p; 0) for a > 0. The field-dependent part 
Jaa (p; a) — JaV (p~, 0) scales in the same way as JaV (p; 0) 
since here A(p;a) = s°A(sp;sa). A similar analysis as 
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FIG. 4: Occupation number as a function of momentum for 
different times. The solid line corresponds to t = 8000, along 
with fits ~ |p| -4 for low and ~ |p| -3 ^ 2 for higher momenta. 

above yields for J d 3 pJa^ (w, p; <p) — in this regime 

a = 3/2. (12) 

Consequently, simulation results may be fitted by pertur- 
bative estimates down to rather low momenta [J]. 

It remains to show that the nonequilibrium evolution 
approaches non-thermal fixed points. For this we fol- 
low the classical-statistical evolution numerically to late 
times. Fig. 2] gives the occupation number n(t, p) as a 
function of three-momentum at different times. We em- 
ploy N — 4 and a weak coupling A = 10~ 5 in order to 
keep quantum corrections small. The latest time shown 
is represented by the solid line. For a region of momenta 
above p > 1 the evolution is very well approximated by 
a power-law behavior n(t, p) ~ | | 3/2 j n a g reemen t with 
([12"]) . As time proceeds this region grows towards higher 
momenta, which is indicated by the arrow in Fig. |4] 

Going to smaller momenta one first observes a transi- 
tion region as expected from the above discussion. The 
infrared behavior is consistent with a power-law that is 
well approximated by a = 4. The evolution in this re- 
gion is very slow, as can be inferred from comparing to 
the t = 3000 line, and getting an even better possible 
agreement between numerics and analytics would be very 
costly in computational time. From the current data the 
a = 5 solution indicated in (jlip is clearly unfavored. 

Scaling solutions govern the nonequilibrium dynam- 
ics only if classical-statistical fluctuations dominate [3] • 
For initial conditions leading to instabilities the statis- 
tical fluctuations grow large. In this case the system 
approaches non-thermal fixed points exhibiting critical 
slowing down. Non-thermal fixed points are unstable 
with respect to quantum corrections, which eventually 
lead to thermalization. The quantum evolution is dia- 
grammatically described by same topologies as displayed 
in Fig. [5]at NLO. However, pairs of propagators f(p)f(q) 



in classical equations ([5]) or JSj) can be associated to 
f(p)f(q) + {^ 2 /^)p(p)p(q) in the quantum theory @. The 
coupling does not drop out and no universal scaling so- 
lutions appear if quantum corrections dominate. 

Non-thermal fixed points can excessively delay ther- 
malization. In the context of early universe reheating 
a most conservative limit requires thermal equilibrium 
at a temperature of order 10 MeV before Big Bang Nu- 
cleosynthesis. Even this might already rule out some 
very weakly coupled inflaton models [4j, however, more 
realistic models have to be considered including quan- 
tum corrections. In the context of heavy-ion collisions, 
QCD plasma instabilities can lead to a regime with qual- 
itatively similar scaling behavior Classical simula- 
tions indicate possible fast isotropization due to instabil- 
ities in about 1-2 fm/c for low momenta of less than 1-2 
GeV |l2j. However, subsequent scaling behavior leads to 
large isotropization times for higher momentum modes 
in this case, which seems incompatible with experimental 
findings. So far, the role of quantum corrections for the 
lower occupied high-momentum modes were not taken 
into account. The question of whether it is possible to 
find also an effective weak coupling in QCD, which may 
facilitate an analytic quantum description, is exciting. 

We are indebted to Szabolcs Borsanyi for many discus- 
sions and computational help. 
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